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This is a semi-annual status report to the National Aeronautics and Space Administration (NASA) 
concerning NASA grant NAG 5-495. This grant was awarded to Dr. C. Megan Urry of the Mas- 
sachusetts Institute of Technology in response to a proposal, entitled Coordinated Observations of X- 
Ray Bright BL Lacertae Objects, to use the International Ultraviolet Explorer (IUE) satellite. The grant 
was awarded on 3/15/85; this report covers the first six months of its duration, or March 15, 1985 
through September 15, 1985. 

In April 1985, I had an IUE observing run at the Goddard Space Flight Center (GSFC) in Green- 
belt, Maryland. Two BL Lac objects, Mrk 180 and Mrk 501, were successfully observed during two 
IUE shifts. Concurrently with the observing run, these and several other IUE images of BL Lac objects 
were analyzed at the IUE Regional Data Analysis Facility (RDAF), which is also at GSFC. The dura- 
tion of the stay was approximately one week. Funds firom grant NAG 5-495 supported this travel. This 
was the only expenditure during the period March 15 to September 1985 as no observing runs were 
scheduled in the remainder of 1985. 

In the same time period, my co-investigators and I worked on a paper with Japanese collaborators 
involving IUE observations of the BL Lacertae object Mrk 421. This paper has now been submitted to 
the Astrophysical Journal, and a copy of the preprint is attached to this report We anticipate that the 
observations of Mrk 180 and Mrk 501, also undertaken as part of collaborative broad-band spectral pro- 
grams, will soon be incorporated into additional publications. 

We have been working with several other IUE users and the staff of the RDAF to develop a new 
program for analyzing IUE spectral images. This method, which extracts spectra from IUE images 
based on Gaussian fits to spectral crosscuts, will be described in a future report, after the software 
development is completed. 

This concludes the first semi-annual report Total direct expenditures for the grant were approxi- 
mately $500 dollars. 
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ABSTRACT 


Simultaneous multifrequency observations of the BL Lac object 

Mkn421 covering radio through X-ray wavelengths were performed on two 

occasions separated by five weeks in 1984 January and March, and each 

observation was coordinated for about a week. We obtained composite 

multifrequency spectra of the central nonthermal component at the two 

epochs after subtracting the optical and infrared light of the 

underlying galaxy. The spectra show the gradual steepening toward high 

frequency; the power law indices are roughly ''O.l, M).6, and 'H.O for 

radio, infrared-optical and UV bands, respectively. The UV and 

optical-infrared fluxes decreased by about 20 percent in five weeks, 

while the radio flux remained stable. The X-ray flux decreased by a 

factor of ^2 and the change was more pronounced at hard X-rays, which 

suggests that X-ray emission possibly consists of two components. The 

degree of polarization at the optical band varied on the time scale of a 

few days while the position angle remained unchanged. 

Physical parameters of Mkn421 are discussed in terms of the 

synchrotron self -Compton model. Taking the spectral turnover between 

infrared and radio for synchrotron self absorption, the radio emission 

originates in a more extended region than the infrared to X-ray 

-2 

emission, the source size of which should be less than about 10 

milliarcseconds . Relativistic beaming is required if the angular size is 

_3 

smaller than a few times 10 milliarcseconds. A possible explanation of 
the spectral change during the two, epochs is also discussed. 
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I. INTRODUCTION 


The central problem of active galactic nuclei (AGN) is the emission 
mechanism of the continuum radiation at radio through gamma-ray 
frequencies, and the ultimate energy source of that radiation. 
Apparently these AGN somehow generate very high luminosities (assuming 
cosmological distances and isotropic emission) inside very small volumes 
(as deduced from variability timescales). BL Lac objects are 
distinguished from other classes of AGN by their rapid variability, high 
degree of polarization, and general lack of discrete spectral features. 
These properties are useful diagnostics of the radiation processes. 
Simplified models such as the synchrotron self-Compton model with or 
without relativistic beaming have been proposed to explain many of the 
observed characteristics of BL Lac objects and other AGN. A high degree 
of polarization in the infrared and optical region may be evidence of a 
synchrotron process, a hypothesis that can be strengthened if correlated 
temporal changes in radio, infrared, and optical bands are detected. 
Many AGN have been observed simultaneously in various frequency bands 
from radio to X-ray frequencies in order to search for correlated 
variations within an individual spectrum (Kondo et al. 1981; Bregman et 
al. 1982, 1984; Worrall et al. 1982, 1984a, 1984b, 1984c; Wills et al. 
1983; Mufson et al. 1984; Glassgold et al. 1983; Hutter and Mufson 
1985). The absence of discrete spectral features, often a nuisance when 
determining the physical state of the emitting region, in this case 
provides a conveniently uncontaminated broad band spectrum. 


Mkn421 was first noted as a blue excess object which turned out to 
be an elliptical galaxy with a bright, point-like nucleus. The redshift 
of the galaxy was measured to be z=0.03 (Ulrich et al. 1975) but the 
spectrum of the nucleus was featureless. It also showed the optical 
polarization and flat radio spectrum, common to BL Lac objects (Maza et 
al. 1978; O'Dell et al. 1978a). The angular size of the radio source 
associated with Mkn421 is less than 1 milliarcsecond (mas) at 6cm 
(Weiler and Johnston 1980; Baath et al. 1981). The ultraviolet emission 
from Mkn421 has been monitored over the past seven years with the 
International Ultraviolet Explorer (IUE) satellite. Nine IUE spectra 
obtained between 1978 and 1981 showed some variability, but no or only 
very weak correlation between flux level and spectral index (Ulrich et 
al. 1984). Since it is one of the brightest BL Lac objects known in the 
X-ray band, Mkn421 has been observed fairly often with various X-ray 
instruments (see summary of observations and references in Urry et al. 
1985). X-ray observation of Mkn421 with HEAO-1 satellite detected a 
change in spectral index from 2.1 to 3.9 and the disappearance of the 
hard component above lOkeV by comparing with earlier observations with 
0S0-8 and SAS-3 satellites (Mushotzky et al. 1979). Because of the lack 
of simultaneity among the observations at radio through X-ray 
frequencies, correlations among intensity changes have been poorly 
studied . 

We made simultaneous multifrequency observations of the BL Lac 
object Mkn421 at radio, infrared, optical, ultraviolet, and X-ray 
frequencies on two occasions in 1984 January and March. In § II we 
summarize observations performed at each frequency band and we describe 
the composite spectrum thus obtained at each of the two epochs. In § III 


interrelations among various observed bands and the comparisons with 
previous observations are discussed. In § IV an interpretation of those 
results in terms of the synchrotron self- Compton model is attempted. 
Conclusions will follow in § V. 

II. OBSERVATIONS 

The observation was conducted at the radio, infrared, optical, 
ultraviolet and X-ray bands on two occasions separated by five weeks, 
1984 January and March. The coverage of the observed frequency bands and 
the instruments utilized in the present observations are summarized in 
Table 1. 


Table 1 


a) X-ray observations 

X-rays were observed with two sets of the scintillation 
proportional counters (SPC-A and SPC-B) on board the TENMA satellite 
(Tanaka et al. 1984; Koyama et al. 1984). The observed energy range was 
1.5 to 35keV, but above about lOkeV, the flux was too faint to 
discriminate from the background. The day-to-day variations of the 
integrated energy flux over the energy range 1.5-10keV are depicted in 
the lowest panel of Figure 1. The statistical accuracy of day-to-day 
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data was not good enough to obtain the energy spectrum. Then we summed 
up the data for each observation period; the time-averaged energy 
spectra observed in January and March are shown in Figure 2. 


Figures 1 and 2 


The flux decreased by a factor of about two from January to March. 

Particular attention has been paid to the spectral analysis, since for 

faint sources such as Mkn421 the spectral shapes depend critically on 

background subtraction. The properties of the background spectrum have 

long been studied carefully using extensive data sets compiled from 

TENMA (Koyama et al. 1984). Before summing the data from SPC-A and 

SPC-B, the spectral analysis was independently carried out for each data 

to check the consistency of the results. The spectra obtained 

independently from SPC-A and SPC-B were combined to improve statistics, 

after confirmation of consistency. The best fit power law index for the 

January and March observations over 1.5-10keV is 1.1±0.2 and 1.8±0.3 

-2 2 

with the normalization factor (1.01±0. 17)xl0 keV/cm sec keV and 
-3 2 

(8.7±2.4)xL0 keV/cm sec keV at IkeV, respectively. 

It is, however, to be noticed that the single power law fit is only 

2 

marginally acceptable, since the reduced x value is 2.3 (degree of 
freedom; 11) and 2.2 (degree of freedom; 7) for the January and March 
data, respectively. The deviations from the single power. law models are 
evident above 4keV. The spectrum obtained in January is flatter above 
4keV as compared to March observations. Below 4keV, on the other hand. 
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the spectral index is 1.6±0.2 and 1.7±0.3, respectively. While the 
intensity varied by a factor of two. 

We notice that the EXOSAT satellite also observed Mkn421 in the 
period 1984 February 1 to 6, one week after our January observation 
(Warwick et al. 1985). The intensity and spectral index are intermediate 
between those of our January and March observations. 

b) UV observations 


Ultraviolet observations were made by the short wavelength prime 
(SWP) and the long wavelength prime (LWP) cameras of the IUE satellite 
on 5 occasions during the present program. The results are summarized in 
Table 2. The integrated fluxes are plotted in the fifth panel of Figure 
1. The spectra were extracted from the line-by-line file using a 
customized narrow slit to maximize the signal-to-noise. Data points 
contaminated by reseaux or cosmic ray hits were removed, and the iata 

c 

were then binned in intervals of approximately 25A. The binned data 
points were weighted according to the rms deviations within the bins and 
were fit to a power law spectrum using a least-squares procedure for 
1200-3000A range. 


Table 2 


It is to be noted that the integrated ultraviolet flux decreased by 
2CK30 percent from January to March and the spectral index a shows a 


? 



slight steepening ( a change of a roughly 0.1). 


c) Optical observations 

Optical observations were made at the Okayama Astrophysical 
Observatory, the Dodaira Station, and the Mt. Lemon Observatory. Five 
telescopes were used with various diaphragm sizes, which are listed in 
Table 3. 


Table 3 


Mkn421 is found at the core of a giant elliptical galaxy and, in 
particular, optical and infrared fluxes are contaminated by galaxy 
emission. Therefore, we must subtract the galaxy component to obtain the 
flux from the central compact core. The distribution of the surface 
brightness of the galaxy was studied by Kinman (1978), Mufson et al. 
(1980), and Hickson et al. (1982) but some discrepancies remained among 
their results. Kikuchi and Mikami (1986) also made the multi-aperture 
photometry of Mkn421 and obtained a result similar to that by Kinman 
(1978). In the present study, the subtraction of the contribution from 
the underlying galaxy has been made using the-results by Kikuchi and 
Mikami (1986). The corrections for the apertures used in the 

t 

observations are also tabulated in Table 3. The colors of the underlying 
galaxy are taken to be constant as U-B*=0.46, B-V=1.05 and V-R«0.90 after 
Kikuchi and Mikami(1985) . These colors are consistent with those of a 
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redshifted elliptical galaxy at z=0.03. For R-I, the mean value of 0.78 
for elliptical galaxies (Johnson 1966) is assumed, since the measurement 
of the underlying galaxy has yet to be made at the I band. 

The observed fluxes and the fluxes of the nonthermal component are 
tabulated in Table 4. The day-to-day variation of the flux of the 
nonthermal component at the V band is shown in the fourth panel of 
Figure 1. The nominal errors associated with the observed fluxes are 
typically 2^3 percent but 1 percent for several cases at Dodaira in 
March. However, the derived flux of the nonthermal component shows a 
scatter greater than the observational errors reaching about 10 percent, 
depending on the aperture and/or telescope. A part of this scatter is 
certainly caused by the incompleteness of the subtraction of the galaxy 
contribution. Although we have no reasonable explanation for the 
remaining scatter, we suspect that the seeing condition and the 
scattered light of a nearby 6 mag star affect the measurement of the 
extended source and of the sky brightness, respectively. 


Table 4 


Continuum spectra for both observing periods are shown in Figure 3 
as well as those in the infrared region. In -the January data of the 
fluxes in the B and V bands, we find a systematic difference between the 
data obtained by the 100cm telescope at Mt. Lemon and those by the 188cm 
telescope at Okayama, although the spectral indices coincided well with 
each other. The difference amounts to 7 percent, and the origin of the 
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difference is unknown. For the spectra in the March epoch, we used the 
results after March 4, since the optical flux abruptly decreased around 
March 3 but stable thereafter. 


Figure 3 


The average V band intensities in January and March were 8.9±0.5mJy 
and 7.8±0.5mJy, respectively. The spectral index in the UVBR region 
varied from 0.65±0.15 to 0.76±0.10 from January to March. Slight 
darkening and spectral steepening are barely significant, if we consider 
only the data in the optical region. It should be also noted that, in 
early February, a brightnening to 10.4±0.3mJy in the V band with a 
spectral index of 0.72±0.06 has been observed at Dodaira. 

Optical polarimetry was also made simultaneously with the 
photometry by the 188cm telescope at Okayama, and by the 91cm telescope 
at Dodaira. The polarization of the nonthermal component showed no 
wavelength dependence in the UBVR region throughout the whole observing 
period, provided that the light from the underlying galaxy was 
unpolarized. The polarization properties of the nonthermal component 
through the U to K bands averaged over March 8-12 are given in Table 5. 
No wavelength dependence of the degree of polarization was found in the 
optical and near infrared regions. 


Table 5 
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The variation of the degree of the polarization is shown in the 
first panel of Figure 1. The degree of polarization showed a 
considerable change during the March epoch; from February 29 to March 5 
the degree of polarization was as high as with a trend of the 
smooth weakening and on March 6 it declined abruptly to 4v5%, while the 
position angle remained constant at 170M80 degrees. This abrupt change 
of the polarization seems to be correlated with the flux decrease in the 
same period . , 


d) Infrared observations 

Infrared observations were made with the 3.8m telescope at United 
Kingdom Infrared Telescope (UKIRT) at Mauna Kea, Hawaii and the lm 
telescope at the Agematsu Infrared Observatory. At Ageiratsu we observed 
Mkn421 on January 25 and 27 and March 7 and 8 with a rectangular 
aperture of 24" in the right ascension by 60" in the declination and 40" 
chopping throw at the J, H and K bands. The calibration was made by 
using £UMa as a standard star and SA0062387 as a reference star. UKIRT 
observations were made on February 4 and 6 at the J, H, K and L bands 
with a 18" aperture using GL299 as a reference star. On March 8-12, the 
J, H, K and L' bands observations were made .with 4.8", 7.8" and 11.8" 
apertures, on March 11-12 the M band observation with a 4.8" aperture, 
and on March 8-10 the N band observation with a 6" aperture. On January 
30, photometry at J, H, K, and L bands was made with a 8" aperture using 
the NASA infrared telescope, at Mauna Kea, Hawaii. 
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The observed data are tabulated in Table 6. The infrared flux 


decreased by L3±4 percent from January to March comparing the Agematsu 
observations at two epochs. As in the optical observations, the 
contribution of the underlying galaxy was subtracted to obtain the 
fluxes of the central nonthermal component for J, H, K, L (L 1 ) bands. 
Since the observation with the largest beam size of 24"*60" does not 
agree with the brightness of the model galaxy, we obtain contrariwise 
the brightness of the underlying galaxy for the 24"x60" beam size by 
subtracting the nonthermal component derived from the observations with 
smaller beam sizes on March 8, and apply it to the photometry in the 
January epoch. The infrared colors of the undelying galaxy are obtained 
as J-H=0.69 and H-K=0.40 from the observations in March, which become 
almost consistent with those of typical elliptical galaxies (Frogel et 
al. 1975; Persson et al. 1979; Glass 1984) taking account of the 
K-corrections (Frogel et al. 1978). We assumed K-L=D.06 (Aaronson 1978) 
and V-K=3.40 (Aaronson 1978; Persson et al 1979), consistent with those 
of typical elliptical galaxies including K-corrections. For the M and N 
bands the contribution of the galaxy is safely neglected. The resulting 
fluxes of the nonthermal component are summarized in the bottom panel of 
Table 6 and the flux at the K band is plotted in the third panel of 
Figure 1. The spectra obtained for respective epoch are shown in Figure 
3 along with the optical spectra. The spectral index in the infrared 
band was 0.5±0.2 and 0.57±0.11 for January and March, respectively. The 
conversion of the magnitude to the absolute flux has been made using 
coefficients tabulated in Table 6, on the basis of the calibration by 
Blackwell et al. (1983). 


74 


Table 6 


The polarization measurements were also made on March 9, 11 and 12 
at the J, H and K bands by UKIRT. Those results are tabulated in Table 5 
together with the optical polarization. The degree of polarization and 
the position angle of the nonthermal component in the infrared are 
4.8±0.4 percent and 180±5 degrees, which are fully consistent with those 
at optical bands. 


e) Radio observations ( 

Radio observations were made with the 45m telescope at the Nobeyama 
Radio Observatory and the 26m telescope at the University of Michigan. 
At Nobeyama, flux measurements at 10, 22 and 43GHz were made, while at 
Michigan flux was measured at 4.8, 8.0 and 14.5GHz. At Nobeyama the 
calibration was done by using 3C286 as the primary calibrator of which 
fluxes are 4.7, 2.6 and 1.7Jy at 10, 22 and 43GHz, respectively, and 
3C274, 4C39.25 and 1144+40 were observed as references. No significant 
daily variations above 10 percent level were found, so that only the 
averaged data over each epoch are used . The results are listed in Table 
7, from which one sees that radio fluxes are rather stable and 560, 520 
and 480mJy at 10, 22 and 43GHz, respectively, throughout January and 


March observations. 


Table 7 


The data obtained at Michigan are consistent with those at 
Nobeyama, 670, 610 and 550mJy at 4.8, 8.0 and 14.5GHz, respectively. In 
the second panel of Figure 1, the data at 14.5GHz are plotted. 

f) Composite spectra 

Figure 4 shows the composite energy spectra emerging from the 
central nonthermal component of Mkn421 in 1984 January and March. These 
composite spectra are -similar to those of other BL Lac objects, i.e., a 
power law spectrum with a gradually increasing spectral index toward 
high frequencies. In the January epoch, the spectral indices are M).l, 
'''0.5, ''-0.7, *vl.O and 'v-l.l for radio, IR, optical, UV and X-ray 
wavelengths, respectively, while in March these a;e ■■‘'0.1, ^. 6 , v 0.8 , 
M . 1 and ^1.8. 


Figure 4 


III. COMPARISON WITH PREVIOUS OBSERVATIONS 

At radio wavelengths, VLBI observations (Baath et al. 1981; Baath 
1984) revealed an unresolved core of 240±40mJy with a size less than 
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0.3mas as. well as a resolved component of 270±50mJy and 3.4mas size at 
5GHz. Other eVLBI observations (Weiler and Johnston 1980; Wehrle et al. 
1984) also -..revealed a compact component of Mmas size and VJOOmJy 
strength both at 5GHz and 2.3GHz. There also exists a 'vO.Uy extended 
component- of a few arcmin size (Kapahi 1979; Ulvestad, Johnston and 
Weiler 1983)6. Our flux measurements at radio frequencies included all 
these components. But the total flux of Mkn421 was rather stable at 
centimeter, wavelength and our value is consistent with previously 
reported vialues (Margon, Jones, and Wardle 1978; O'Dell et al. 1978a). 
Therefore the flux of the core component is probably about 300 mJy. 

O'Dell et al. (1978a) showed that the. radio spectrum extends to at 

least 90GHz Vwith a flat shape. The extrapolation of the flat radio 

spectrum of c 300mJy intersects the extrapolation of the IR-optical 

12 

spectrum at ^around 10 Hz . 

—The .optical fluxes were faint compared to the previous observations 
(Hagen-Thoifioret al. 1983). Also the infrared fluxes of the nonthermal 
component -were in a faint phase among previous observations at the K 
band (UlrichLtet al 1975; Allen 1976; Joyce and Simon 1976; O'Dell et al. 
1978b; Bailyiiet al. 1981; Impey 1983). Thus we conclude that Mkn421 was 
in a dark phq.se on two occasions of our coordinated observations. 

The position angle of the optical polarization in the present 
observation agrees with the dominant direction in the long term behavior 
(Hagen-Thorn ret al. 1983). The polarization of -Mkn421. observed by Baily 
et al. (1981)6 in 1980 December is nearly equal to that in the present 
observatioitdin the position angle and the degree of polarization, though 
Mkn421 was in :a brighter phase then. However Bailey et al. claimed that 
the polarizaitlon degree decreases as the wavelength increaces, while our 
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degree of polarization shows no wavelength dependence. We think it is 
not due to observational errors but due to the difference in the models 
of the underlying galaxy. Their model galaxy has the brightness by 
0.3mag smaller than that we used in the V band, which brings a 
decreasing degree of polarization as the wavelength increases through 
the dilution effect. Therefore, a decreasing polarization of Mkn421 with 
wavelength is questionable. 

It is to be noted that there exists a spectral break between 
IR-optical and UV bands at around lO^Hz. This spectral break or 
steepening seems to be common to BL Lac objects. 

The UV light curve for the period 1978-1981 was analyzed by Ulrich 
et al. (1984), and they have shown that the flux varied on the time 
scale of a few years. In 1979 November the UV flux was in a minimum 
state with the spectral index of 1.0. Both the flux and the spectral 
index are similar to those of the present observation in 1984. Since the 
flux steadily increased and the spectrum hardened from 1979 to 1981, 
Mkn421 might have turned to a decreasing and softening phase in 
1982-1983. 

The X-ray flux of this observation is significantly lower than 

those of previous observations by SAS 3 (Hearn et al. 1979) and HEAO 1 

(Mushotzky et al. 1979) by an order of magnitude. At those epochs X-ray 

fluxes were higher than or on the extrapolation from the UV spectrum. 

The extrapolation of the UV spectrum toward shorter wavelengths should 

be compared to the X-ray spectra obtained by TENMA. The extrapolation 
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from the UV spectrum lies at 5.47 and 2.26vJy at 10 Hz at January and 
March epochs, respectively, while the observed fluxes were 1.1 and 
0.51VJy, respectively. Then the ob.served X-ray fluxes are lower than the 



extrapolation from the UV spectrum by an order of magnitude in the 

* 

present observation. The spectral slope between UV and X-ray is 1.3M.4 
and apparently steeper than the UV slope of M. On the other hand, the 
X-ray spectral index changes from M at the January epoch to ^2 at the 
March epoch. Below 4keV the X-ray spectral index was about 1.5 for both 
periods . 

This result implies that X-ray emission is represented by two 
components. These variations of flux and spectral shape at five week 
interval are very important feature. HEAO-1 observations in 1978 May 
showed a remarkable spectral change from 0S0-8 observations in 1977 May 
(Mushotzky et al 1979) . It has been suggested that the X-ray spectra of 
BL Lac objects have two components, a soft and a hard component (Riegler 
et al. 1979; Worrall et al. 1981). We may infer that the hard component 
varies more violently than the soft component. 


IV. IMPLICATIONS ON THE MODEL 
a) Source parameters 

The continuum emission from BL Lac objects, may be due to 
synchrotron radiation and/or inverse Compton scattering (Blandford 1984; 
Konigl 1981). In this subsection we derive physical parameters of the 
source using the data obtained in the present observations within the 
context of the homogeneous and spherical symmetric synchrotron 
self-Compton (SSC) model (Gould 1979). The SSC model assumes the 
presence of relativistic electrons < with power law distribution, N(y)*Y ^ 
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(Y is the Lorentz factor and p is the power law index), which emit 
synchrotron radiation having the power law spectrum ^ S y n “ v ^ Since 
the synchrotron absorption coefficient is proportional to v t ^ e 
self absorption turnover is expected at a certain frequency v (v GHz) 
with the corresponding flux density S^Jy. 

Common to those BL Lac objects for which the optical and UV fluxes 
were observed, a break or bending of the spectrum appears at the break 
frequency v MO^Hz. The acceleration rate of relativistic electrons 
responsible for synchrotron radiation above may be smaller than the 
energy loss rate, making the energy spectrum steeper. Here the Lorentz 
factors of electrons which emit the synchrotron radiation at the 
turnover and break frequencies are denoted by Y and Y, , respectively. 

3 D 

Soft x-rays in January epoch and the whole X-ray in March epoch may be 
produced by the synchrotron radiation by electrons in the high energy 
tail. The relatively flat hard X-ray component as seen in the January 
spectrum (see Figure 2) may be attributed to the self-Compton photons, 
i.e., inverse Compton scattering of synchrotron photons off relativistic 
electrons. This self-Compton flux has a power law spectrum with the 
same spectral index as the synchrotron spectrum with 


S « v " ( P" 1)/2 e " 2(p+2) v _(3 P +7)/2 S (p+3) (6/(l+z))“ (p+3) ln(v /v ), (1) 

com X * x * b a 


where is the frequency of a scattered photon, 6 is the kinematical 
Doppler factor in the relativistic beaming model, z is the cosmological 
redshift, and the apparent angular diameter of the source is O^mas. 

Eenergy densities of magnetic fields u ma g and synchrotron photons 
u , the lower limit to the enejrgy density of relativistic electrons 
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Tn< rv 

u . , and y are given as 
rel a 

U mag* KWX'^HI+z)) 2 , 

"syn“ e *" 2 ' , * s *(' s /( 1+z )>’ 4 - 

min n -9 -7_ 4 /x//1 , >.-5,. , .2, -1 

U rel“ ®* V * S * (6/(1+z)) (1+z > d L * • 

and 

V 0*‘ 2 v*" 2 S*(6/(1+ Z ))' 1 , 


( 2 ) 

(3) 

(4) 

(5) 


respectively (Marscher et al. 1979; Marscher 1983). Here d^Gpc is the 
luminosity distance. The proportional coefficients in Equations (1)'\'(5) 
are complicated functions of p. 

Since the X-ray flux by: the self-Compton photons should be lower 
than the observed X-ray flux, we obtain the lower limit to the beaming 
factor 6 from Equation (1) as. 


lnCv^/v^) l/(p+3) 

6 min“ (1+Z)S * ( (p-1) /2 2 (p+2) (3p+7)/2 } 

X * * 


( 6 ) 


The synchrotron self absorption turnover was not directly observed, 

3 

but presumably v^. lies at around MO ,- - the intersection of 
extrapolations from radio and IR spectra. Assuming that the flux at the 
turnover is not much different from that of the smallest unresolved 
radio component (0.3Jy with an angular size less than 0.3raas), we adopt 
S^=0.3. The turnover frequency .is determined so that the synchrotron 


21 


flux gives the observed UV flux at lO^Hz (6.5mJy for the January epoch 

*■ 

and 5.1mJy for the March epoch) with the spectral index ot=0.625 (p=2.25) 

the well fitted value to the IR-optical observations. The source 

parameters are calculated as a function, of 0^, which are shown in Figure 

5 for the January data. We notice that there are two critical values of 

0*, referred to as 0^ and 02* For 0^>0^, Y a is less than unity. Since 

the synchrotron flux with a power law spectrum is valid for y>>l, thus, 

0^ provides an optimistic upper limit to the angular diameter. While for 

0.<0„, the expected inverse Compton flux at v deduced from the observed 
x z X 

synchrotron flux exceeds the observed flux if 6=1, i.e., no reletivistic 

bulk motion of the source is invoked. Conversely such relativistic 

effects are required if 0* is found virtually smaller than In the 

first and second columns of Table 8, the results for the January and 

18 

March data are tabulated for 6=1 and 0 A =02» taking v x = 10 Hz. 


Figure 5 and Table 8 


Ulrich et al. (1975) reported the absorption redshift of the 
underlying galaxy, z=0.03, implying the distance of 90h *Mpc; lmas 
corresponds to 0.44h ^pc and the light crossing time of the source is 
t^ c =4.5 x 10^8 Jk h ^sec, where h is the Hubble parameter in' units of 100 km 
s * Mpc ^ . In Table 8 the cooling times, t and t, , of electrons with 

3 u 

Lorentz factors of Y fl and Y^ • are given respectively, as well as t^. 
The cooling time of electrons responsible for the IR-optical radiation 
is shorter than the light crossing time, implying in situ acceleration 



of the relativistic electrons. Also from Table 8 one sees that the 

energy density of magnetic fields, u , dominates over that of 

mag 

relativistic electrons, noting that is not so different from true 

energy density of relativistic electrons since y is relatively small. 

a 

This property is consistent with the observed behavior of polarization 
in the IR-optical bands that the position angle remained constant 
although the degree of polarization and flux varied with time. 

Since VLBI experiments have shown that 'ulmas structure exists at 
5GHz, radio structure should be more extended than IR-optical-UV 
component. This is consistent with the model of flat spectrum radio 

sources which can be expressed in terms of inhomogeneous synchrotron 
source (Marscher 1977; Konigl 1981). In this model the apparent angular 

size is inversely proportional to the frequency in the flat portion of 

12 -3 

the spectrum, and the size at 10 Hz becomes ^5*10 mas, consistent with 
above analysis. 

Mkn421 was rather faint in the epoch of the present observation and 

44 -1 

the total luminosity was about 10 erg s . At the most luminous epoch 

in the past it was about 10 times brighter. If we adopt that the 

ultimate energy release occurs as a result of accretion onto the 

45 -1 

supermassive black hole and 10 erg s corresponds to the Eddington 

luminosity, the mass of the black hole is 10^M@ and the accretion energy 

13 -5 

is reprocessed in the size of 10 cm which amounts to $*=10 . Thus the 

2 

size of the synchrotron source is about 10 times larger than the size 
of the ultimate energy release region. 

(2) Spectral variation 
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' 




Our January and March observations of Mkn421 suggest that the 
self-Compton flux corresponding to the relatively flat hard ' X-ray taii 
seen in the January observations decreased by more than factor ^2 during 
five weeks, while the synchrotron flux decreased by 1 5^2 5 percent in the 
IR, optical and UV bands in the same interval. The radio flux remained 
essentially constant during those two observational periods. It is , 
however, rather difficult to explain the change of the spectra obtained 
in those two epochs with the homogeneous and spherical symmetric SSC 
model as will be discussed below. 

As seen in Equation (2)-(5), the physical source parameters are 

determined as functions of p, 9*, S^, v A and 5. Comparing the composite 

spectra of the two epochs we may approximate that and p remained 

constant, and that decreased by 35 percent resulting in a decrease of 

the synchrotron flux by 22 percent. Let us first consider that alone 

decreased, and 0^ and 6 remained constant. Physically the decrease of 

is a consequence of the decrease of u and the increase of u 1 ’"'?. In 

mag rel 

this case the self-Compton flux increases by a factor of 14 since 

S c v. (^p+7) f r0 m Equation (1), on the contrary to the observed 
com * 

spectral change. 

Next let us assume that in addition to the decrease in v A the 

“8 

source expands as where 6 is an arbitrary constant. From 

8 58-6 875 

Equations (1) ^ com ccv * * * for P = 2.25. In order to reduce the 

1 8 

self-Compton flux at 10 Hz by a factor of 2.2, from l.lyJy in January 


to 0.51pJy in March, we have 3>1.04. In this case, however, 


. 10-86 . win. 96-7 . ^ 

u : change as u <*v. and u ,«v, . As a result, 

rel mag * rel * 


magnetic energy increases, since it chages as 6* u 


10-136 


u and 
mag 

the total 
. Thus, in 


order to explain the observed spectral variation, the increase in the 


magnetic field energy is necessary despite the expansion of the source. 

In the case of adiabatic expansion, u is proportional to 0. \ 

mag * 

0 21 

implying 8=5/6, which results in S «v. * . Therefore, the decrease of 

com * 

the self-Compton flux is less than that of the synchrotron flux in 

contrast to the observational fact. 

Finally, we consider the change of the relativistic beaming factor 

-4 

along with the adiabatic expansion of the source. For u “0^ , 

10 2 -12 

Equation (2) implies 6 *9^ , from which we obtain 

S «v. ^ In orc j er to decrease the self-Compton flux by 

a factor of 2.2, the kinematical Doppler factor is required to increase 

by 20 percent within five weeks and the source would expand during the 

above time interval by 34 percent. Accompanying this and U rel wou ^ 

decrease by a factor of 3.2 and 2.3, respectively. Thus the increase in 

the relativistic motion of the source can account for the observed 

property of the spectral variation, separated by five weeks, although 

some acceleration mechanism must be. invoked. In the third column of 

Table 8, the source parameters in the March epoch are shown assuming 

that the source in January was in a state of an< ^ $ = 1* 

Since the hard X-ray spectrum was very steep in March, we obtain 

more severe limit if we apply the same argument at higher v . In the 

A 

fourth column of Figure 8 are tabulated the numerical results when we 
1 8 

take v^=1.4*10 Hz. 

We have not applied the inhomogeneous jet model presented by Konigl 
(1981), since some artifacts are invoked in this model; the magnetic 
field and relativistic electrons are distributed so as to determine fit 
the whole range of the composite spectrum. However, the analyses of our 
observations with such models are also awaited. 


V. SUMMARY 


The results of the present observations are summarized as follows. 

(1) Mkn421 was in a relatively faint state in 1984 January and March at 
frequencies from IR through X-rays. 

(2) From January to March, Mkn421 decreased its flux in the IR to X-ray 
bands. X-ray emission decreased by a factor of 2VJ while IR-UV radiation 
decreased by 15^25 percent. Radio emission remained essentially stable. 

(3) The change in the X-ray spectral shape suggests that X-ray emission 
has two components; the soft component is probably the synchrotron 
radiation and the hard component was more violently variable and may be 
ascribed to other mechanisms such as the inverse Compton scattering of 
synchrotron photons. 

(4) The break was observed at a frequency of around lO^Hz; The power 
law index of IR-optical part is O^^.S, while for the UV range it is 
M.O. 

(5) At the optical band, the degree of polarization showed a factor 2 
variation on the time scales of a few days, while the position angle 
remained constant during the present observations. The degree of 
polarization in the IR-optical region shows no wavelength dependence. 

(6) The radio spectrum is flat and the extrapolation from the radio 
spectrum intersects the extrapolation from IR at a frequency of a few 
times lO^Hz, below which the synchrotron self absorption plays an 
important role. 

(7) Applying the homogeneous synchrotron source theory to the observed 

results, physical source parameters are deduced. Source size should be 

_2 

less than about 10 mas. If the source size is less than a few times 


Z(> 


_3 

10 mss, relativistic beaming is required. The source is magnetic field 
dominated and the relativistic electrons with relatively low energy are 
responsible for the synchrotron radiation. The cooling time of electrons 
is shorter than the light crossing time, implying that there must be in 
situ acceleration of the electrons. 

(8) The spectral variations in the five week interval are hard to 
interptete. One possibility is that the source expanded along with an 
increase of relativistic beaming factor. 
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Table 1 


Organization of observations 


Observatory 

Band 

Observers 

TENMA 

X-Ray 

Mak ino , Tanaka , Mat suoka , Koyama , Inoue , 


1.5-30keV 

Makishima ,Hoshi ,Hayakawa 

IUE 

UV 

Kondo ,Urry,Muf son, Hackney .Hackney 


1200-3000A 


Mt .Lemon 

optical 

U,B,V,R,I 

Wisniewski 

Okayama 

optical 

Kikuchi .Mikami 

and Dodaira 

U,B,V,R 



with polarimetry 

Agematsu 

IR 

J,H,K 

Hiromoto,Nishida 

UKIRT 

IR 

J,H,K,L,L' ,M,N 

Burnell, Brand .Williams, Smith 


with J,H,K polarimetry 

NRO 

Radio 

Takahara , Inoue .Tsuboi .Tabara ,Kato 


10, 22, 43GHz 


Michigan 

Radio 

Aller.Aller 


4. 8, 8.0, 14.5GHz 


3i 



Table 2 


Results of UV observations by IUE 


Date 

Image 

A +) 

(mJy) 

a”> 

x 2 /dof 

Integrated 

*) 

* flux ’ 

Jan. 23 

LWP2669 

6.65±0. 12 

1 . 04±0 .04 

0.71 

6.0110. 11 


LWP2700 






SWP22082 





Jan. 28 

LWP2711 

6.35+0.20 

1.01±0.07 

0.55 

5.8210.11 


LWP2712 






SWP22L28 





Feb. 2 

LWP2732 

6. 77±0. 15 

1 . 13±0. 20 

0.11 

5.8810.18 

Mar . 3 

LWP2881 

5. 24±0. 18 

1 . 10±0.08 

1.09 

4. 6110. 16 


LWP2882 






SWP22398 





Mar. 9 

LWP2915 

5.04±0.21 

1.1210.08 

0.47 

4.3810.14 


SWP22445 






+) the flux density at 3000A 

if) the spectral index for 1200-3000A 

o -11 _2 -1 

*) integrated over 1200-3000A and in units of 10 erg cm s 


3 2s 






Table 3 


List of telescopes for optical observations 

and the aperture correction for the underlying galaxy 


Observatory Diameter Aperture Date 
(cm) (arcsec) 


Aperture Correction 
V magnitude 


Okayama 

91 

37.1 

Jan23 

13.82 



27.9 

Jan29 ,Mar3 ,4 ,6 

14.05 


188 

8.9 

Jan29 ,31 

14.98 



11.8 

Jan28 

14.72 

Dodaira 

91 

18 

Feb7 , 29 

14.39 



13 

Feb. 6, 9, Marl-12 

14.65 

Mt .Lemon 

100 

29.7 

Jan24-Febl ,Mar9 

14.00 


150 

12.4 

Marl-8 

14.68 


Table 4 

Observed optical fluxes and derived fluxes of non-thennal component 


Observed non-thermal 

Date Aperture F(U) F (B) F(V) F(R) F(I) F(U) F(B) F(V) F (R) F(I) 


1984 (UT) 

arcsec 


mJy 




mJy 








Okayama 






Jan23 . 80 

37.1 

8.96 

12.25 

18.69 


7.61 

7.51 

7.22 



25.74 

27.9 

8.40 

11.70 

16.89 


7.31 

7.86 

7.61 



28.80 

11.8 

7.94 

9.91 

13.79 

19.52 

7.36 

7.84 

8.79 

10.60 


29.64 

8.9 

7.73 

9.56 

13.05 

17.81 

7.27 

7.93 

9.11 

10.78 


31.65 


8.40 

10.38 

13.67 

17.48 

7.93 

8.75 

9.73 

10.46 


Mar 3.53 

27.9 

7.38 

11.07 

17.53 


6.29 

7.23 

8.24 



4.82 


7.05 

10.58 

18.02 


5.96 

6.74 

8.73 



6.74 


7.31 

10.29 

16.89 


6.23 

6.45 

7.61 







Mt .Lemon 






Jan24 . 43 

29.7 


11.38 

18.02 



7.36 

8.30 



25.40 



11.07 

17.53 



7.05 

7.81 



26.43 



11.70 

18.52 



7.68 

8.80 



28.33 



10.77 

17.85 



6.75 

8.13 



29.39 



11.38 

17.37 



7.36 

7.65 



31.44 



11.38 

18.02 



7.36 

8.30 



Feb 1.34 



11.70 

18.18 



7.68 

8.46 



Mar 1.35 

12.4 

7.25 

9.73 

14.31 

20.44 27.64 

6.64 

7.58 

9.12 

11.18 

12.31 

2.26 


6.79 

9.38 

13.42 

20.26 27.39 

6.18 

7.23 

8.22 

11.00 

12.06 

3.34 


6.67 

9.04 

13.54 

19.52 27.14 

6.06 

6.89 

8.35 

10.26 

11.81 

4.31 


6.67 

8.88 

13.54 

19.34 26.64 

6.06 

6.73 

8.35 

10.09 

11.31 

5.36 



9.04 

14.05 



6.89 

8.85 



6.30 


6.25 

8.64 

12.58 

18.82 25.45 

5.64 

6.49 

7.38 

9.56 

10.11 

7.26 


6.73 

8.88 

13.67 

19.70 25.68 

6.12 

6.73 

8.47 

10.45 

10.35 

8.30 


6.25 

8.56 

12.93 

18.65 25.92 

5.64 

6.41 

7.74 

9.39 

10.59 

9.29 

29.7 


9.38 

15.27 



5.36 

5.54 




Feb 6.79 

13.0 

8.27 

10.67 

15.23 

Dodaira 

22.08 

7.65 

8.46 

9.89 

12.56 

7.83 

18.0 

9.42 

11.95 

17.55 

25.25 

8.63 

9.14 

10.76 

13.15 

9.78 

13.0 

8.58 

10.86 

15.91 

22.66 

7.96 

8.65 

10.56 

13.14 

29.76 

18.0 

7.72 

10.29 

15.56 

22.43 

6.92 

7.48 

8.77 

10.34 

Mar 1 . 65 

13.0 

7.29 

9.68 

14.22 

20.46 

6 . 66 

7.47 

8.87 

10.94 

2.60 


7.03 

9.26 

13.62 

19.61 

6.40 

7.05 

8.28 

10.09 

4.63 


6.25 

8.40 

12.82 

18.80 

5.62 

6.19 

. 7.48 

9.28 

5.64 


6.36 

8.55 

12.93 

18.80 

5.74 

6.34 

7.59 

9.28 

6.71 


6.30 

8.55 

12.93 

18.80 

0.68 

6 .34 

7.59 

9.28 

7.68 


6.36 

8.55 

12.93 

18.64 

5.74 

6.34 

7.59 

9.13 

8.61 


6.66 

8.86 

13.27 

19.12 

6.03 

6.65 

7.93 

9.60 

10.68 


6.54 

8.78 

13.39 

19.44 

5.91 

6.57 

8.04 

9.93 

11.60 


6.78 

9.02 

13.39 

19.28 

' 6.15 

6.81 

8.04 

9.76 

12.71 


6.48 

8.70 

12.82 

18.64 

5.85 

6.49 

7.48 

9.13 




Table 5 


The polarization of the non-therraal component 
observed in 1984 March 8-12 



U 

B 

V 

R 

J 

H 

K. 

arcsec 

13 

13 

13 

13 

7.8 

7.8 

7.8 

% 

4.0±0.3 

3.3±0.3 

3.010.2 

2.710.2 

2.510.4 

2.910.3 

2.810.6 

% 

4. 3+0. 3 

4.310.4 

4.710.4 

5.0+0.4 

4.410.7 

5.010.4 

4.5H.0 

degree 

177+2 

175+3 

17412 

17513 

17919 

18015 

184+13 


In each column the first, second, third and fourth rows show the 
beam size, the observed degree of polarization, the degree of 
polarization of the nonthermal component, and the position angle, 
respectively. 


if 



Table 6 


Summary of infrared photometric observations 


Date 

Aperture 

J 

H K 

L 

L* 

M 

N 

(UT) 

(arcsec) 

240THz 

182THz 136THz 

86.9THz 

78.5THz 

62.7THz 

28.6THz 



1.25um 

1.65pm 2.2pm 

3.45pm 

3.82um 

4. 78Pm 

10.5um 



1540Jy 

980Jy 670Jy 

290Jy 

250Jy 

165Jy 

40Jy 


Jan25/27 

24 x 60 

46.94 

57.44 

56.24 





Jan30 

8 

25.56 

31.28 

33.89 

32.69 




Feb4/6 

18 

38.7 

49.1 

49.4 

40.4 




Mar7 /8 

24*60 

44.83 

51.91 

54.21 





Mar8 

4.8 

22.06 

26.02 

30.62 


30.90 




7.8 

26.76 

31.86 

35.49 


31.76 




11.8 

31.88 

37.95 

41.50 






6 







58±18 

Mar9 

4.8 

23.52 

26.74 

33.27 


35.15 




7.8 

26.27 

31.28 

34.52 


32.65 




6 







58±12 

MarlO 

7.8 

26.03 

30.71 

32.97 


29.8 




6 







36.5±7 

Marl 1 

4.8 

22.26 

26.26 

29.52 


30.9 

38±6 



7.8 

26.03 

31.57 

34.20 


30.3 



Mar 12 

4.8 

21.65 

26.26 

25.95 


30.0 

29±6 



7.8 

26.27 

31.86 

34.52 


35.5 




Jan25/27 

17.7 

23.9 

24.2 



Feb4/6 

16.6 

22.5 

23.4 28.5 



March8/12 

15.2 

18.1 

21.6 

27.7 33.2 

50.7 


The upper panel gives the central frequency, central wavelength and the 
absolute flux at zero magnitude. The second panel gives the observational 
results and the third panel shows . the fluxes of the non-thermal component in 
units of raJy. The errors are less than 1 percent unless otherwise indicated. 


Table 7 


Summary of Radio Observations 


(a) NRO 45m 

Epoch 10GHz 22GHz 43GHz 


Jan. 0.59±0.03 0.47*^ 0.49+0.03 

March 0.55±0.03 0.55±0.05 0.55±0.1 


(b) Univ. Michigan 26m 


Epoch 4.8GHz 8.0GHz 14.5GHz 


Jan. 0.657±0.034 0.604±0.013 0.540±0.017 

March 0.659±0.036 0.62 ±0.03 0.571±0.020 


The flux density is shown in units of Jy and the errors are the 
standard deviation from the mean of daily flux. 

*) the data are available only for one day-. and the error cannot be 
estimated . 


*7 






Table 8 


The source characteristics deduced from 
the synchrotron self Compton model 


January March 


v* (GHz) 

2.2*10 3 

1 . 5*10 3 

1.5*10 3 

1 .5 X 10 3 

0^ (mas) 

1.6*10“ 2 

2 . 3*10~ 2 

2.3*10 -2 

2.3*10~ 2 

02 (mas) 

1.8*10“ 3 

2.7*10" 3 

2.7*10 -3 

3 . 1*1 0 -3 


6 

1 

1 

1.2 

1.7 

0 A (mas) 

1.8*10” 3 

2. 7*10 -3 

2.4*10“ 3 

2 . 2*10 -3 

t j c (sec) 

7.6*10 4 

1.2*10 5 

1.0*10 5 

9.7*10 4 

Y a 

80 

76 

80 

63 

B (G) 

1.3*10 2 

9.5*10 1 

7.0*10* 

7.9*10 1 

u (ere cm 3 ) 

6.4*10 2 

3.6*10 3 

2.0*10 2 

2.5><10 2 

mag 

u (ere cm ) 

1.2*10^ 

4.0*10 

2.1*10 

5.7*10 

syn 

U rel (erg Cm_3) 

3.1 

1.2 

1.4 

3.9*10'* 

t (sec) 

3.3*10 2 

1.1*10 3 

1.8*10 3 

8.5*10 3 

t (sec) 

1.6*10 L 

3.9*10* 

8 .-4*10* 

3.9*10 2 




Figure Captions 


Figure 1: The light curves of the nontherraal component of 

Mkn421. The first panel shows that of the degree of polarization 
at the V band. The following panels show that of fluxes in the 
radio (14.5GHz), IR (K band), optical (V band), UV and X-ray 
frequencies, respectively. 

Figure 2: The X-ray spectra obtained by the TENMA satellite. The data 
were integrated over each observation period . 

Figure 3: The IR-optical spectra of Mkn421 in 1984 January-February 
(Figure 3a) and in March (Figure 3b) . The observed fluxes are 
indicated by filled circles with straight lines: the beam sizes 
are also shown numerically in the figure. The derived 
nonthermal fluxes are indicated by open circles. 

Figure 4: The composite spectra of of the non-thermal component of 
Mkn421 in 1984 January and March. 

Figure 5: Source parameters obtained by the homogeneous synchrotron 

theory using the January data. The energy density of magnetic 

fields u , that of synchrotron photons u , and the lower limit 
mag J r syn 

to the energy density of relativistic electrons are shown as a 

function of the angular diameter of the source Q^mas. The minimum 


V 


beaming factor 6 , and the Lorentz factor of electrons Y which 
min a 

emit the synchrotron radiation at the self absorption turnover 
frequency are also depicted. 


an 
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